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ABSTRACT 

H and He features in photospheric spectra have seldom been used to infer quantitatively the 
properties of Type lib, lb and Ic supernovae (SNe lib, lb and Ic) and their progenitor stars. 
Most radiative transfer models ignored NLTE effects, which are extremely strong especially 
in the He-dominated zones. In this paper, a comprehensive set of model atmospheres for 
low-mass SNe Ilb/Ib/Ic is presented. Long-standing questions such as how much He can be 
contained in SNe Ic, where He lines are not seen, can thus be addressed. The state of H and 
He is computed in full NLTE, including the effect of heating by fast electrons. The models 
are constructed to represent iso-energetic explosions of the same stellar core with differently 
massive H/He envelopes on top. The synthetic spectra suggest that O.O6-O.14M0 of He and 
even smaller amounts of H suffice for optical lines to be present, unless ejecta asymmetries 
play a major role. This strongly supports the conjecture that low-mass SNe Ic originate from 
binaries where progenitor mass loss can be extremely efficient. 
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1 INTRODUCTION 

Stripped core-collapse supernovae of types lib, lb and Ic (SNe lib, 
lb and Ic) are increasingly well observed, and long-standing is- 
sues about the explosion physi cs dMezzacappal 120051; I Janka et alj 
20071) and the pr ogenitor stars dGeorgy et alj|2009t ISmarttll2009t 



Yoon et al.|[201oh are progressively being resolved. However, one 



possibility for analyses of such SNe has seldom been exploited: ra- 
diative transfer modelling of the H and He features in the spectra 
jUtrobin|[l99o1 : iMaurer et alj|20ld) . Thus, H and He abundances 
have rarely been quantified. The amount of H and He that can be 
"hidden" in Type Ib/c SNe, which do not show clear H or H e lines, 
respectively, is not accurately known (e.g. lSauer et alj|200of) . 

There is ample motivation for studying the H and He layers 
of stripped-envelope SNe. The analysis of abundances and den- 
sities will allow us to validate progenitor models for SNe Ib/Ic 
and address the question of envelope stripping. Information on the 
state of the outer ejecta is also essential to understand early-time 
light curves. At the same time, these layers may play an impor- 
tant role in the SN-gamma-ray-burst (SN-GRB) co nnection: thick 
He or even H shells are thought to quench jets (e.g. iMazzali et al.1 



2008 ) , from which long GRBs supposedly o riginate dWooslevl 



19931 : iMacFadven & W ooslev 1999; Piran 2005) 



Spectral modelling of the He-rich layers of SNe has an added 
complication: the calculation of the state of the He-rich plasma. 
The excitation/ionisation state of the He component must devi- 
ate significantly from a local thermodynamic equilibrium (LTE) 
configuration in orde r for He I lines to be present in SNe lb 
dHarkness et al.|[l987h . Processes leading to large departures from 
LTE have to be explicitly included in the si mulations. The rele- 
vant excitation mechanisms were discussed by |ChugaildT987bh and 
iGrahaml dl988l) . They recognized that He is strongly affected by 
collisions with fast electrons, which result from Compton processes 
with 7-rays from the decay of 56 Ni and 56 Co. 

iLucvl { fl99ll) devised a treatment of excitation and ionisation 
by nonthermal electrons with typical energies >10keV and syn- 
thesised pure He I line spectra on a thermal continuum, under en- 
velope conditions typical for core-collapse SNe. For modelling the 
non-LTE (NLTE) state of He, he devised a treatment of the exci- 
tation and ionisation by nonthermal electrons with typical energies 
>10 keV. Then, he simultaneously calculated the radiative transfer 
and solved the NLTE rate equ atio ns (equations of s t atistic al equi- 
librium). Later. lUtrobir] d 19961) and lMazzali & Lucvl dl998l) treated 
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SNe with He with detailed radiative transfer codes, including NLTE 
calculations. Here, we adopt a similar treatment for He and H and 
present a state-of-the-art code module (subsequently called "NLTE 
solver" or "NLTE module") th at integrates wi th th e radiative tra ns- 
fer / spectral synthesis code of lMazzalU J2000t> and lLucvl dl999l) . 

We compute and analyse spectral models for low- to 
moderate-mass SNe Ilb/Ib/Ic in the photospheric phase, focussing 
on the question of how much H or He can be hidden in a SN with- 
out producing spectral lines. First, we use extensive observations of 
a SN lib (SN 2008ax) and a SN Ic (SN 19941) in order to construct 
envelope models, respectively, by fitting the observed spectra us- 
ing our radiat ive transfer code w ith the NLTE module [abundance 
tomography - IStehle et al .1 J2005h l - In a second step, we use the in- 
ferred abundance and density stratification for the two SNe to set up 
a sequence of SN Ilb/Ib/Ic transition models. The transition models 
are characterised by decreasing H/He envelope masses until only C 
and O remain. We can thus assess the H and He mass required to 
identify these elements unambiguously in the observations. 

We first outline the spectral synthesis code, the NLTE mod- 
ule (with model atoms), and the integration of the two (Sec. [2)- 
We then present the SN lib and Ic models (Sec. [3) and the transi- 
tion sequence, and give limits for the H and He masses (Sec. [4^. 
We further discuss the results in Sec. [5] Finally, we summarise our 
findings and draw conclusions in the context of current progenitor 
models for stripped-envelope core-collapse SNe (Sec. [6). 

2 SPECTRAL MODELS WITH H AND HE IN NLTE 

We use a Monte-Carlo (MC) radiative transfer code to calculate 
mode l spect ra of SNe for our stud y. Our co de has been presented b y 
Lucvl ill 9991) and iMazzalil d2000h [see also lAbbott & Lucvl d 19851) . 
Mazzali & Lucvl dl993l) l. and incl udes the possibility to specify a 
multi-zone chemical composition dStehle et alj |2005). Apart from 
the abundances, the input parameters are a density distribution, and 
- for each epoch at which we calculate spectra - the emergent lu- 
minosity Lboi and the "photospheric velocity" u p iQ. 

2.1 Original code, integration of the NLTE module 

2.1.1 MC transport and MC estimators 

For a given epoch, the MC code computes the radiative transfer 
through the SN ejecta above a boundary ("photosphere"), below 
which the ejecta are considered to be optically thick. From the pho- 
tosphere (at « p h), thermal radiation [I„ = B v (T p h)] is assumed to 
be emitted into the atmosphere. The radiation is simulated in the 
form of energy packets, which undergo Thomson scattering as well 
as line excitation-deexcitation processes. Line transfer is treated in 
the Sobolev approximation. After an energy packet has excited an 
atom, it can be re-emitted as a whole in a different downwards tran- 
sition (photon branching). The downwards transition is chosen ran- 
domly, ensuring correct sampling of the emission in the limit of 
many excitations. This "indivisible packet approach" enforces ra- 
diative equilibrium. The state of the atmosphere is assumed to be 
stationary, as photon transport through the simulated region is much 
faster than changes in the SN luminosity or in the ejecta radius. 
In the course of a MC run, values describing the radiation 

1 In the context of homologous expansion, which applies at the epochs 
considered here, radius r and velocity v can be used interchangeably as 
coordinates (r = v ■ t, where t is the time from the onset of the explosion). 



field [MC estimators - iMazzali & Lucvl dl993h : lLucvl dl999l) l are 
recorded at the midpoints of the shells into which the envelope is 
discretiserfl For the present work, two estimators are particularly 
relevant. One is the angle-ave raged intensity in the blue wing of 
each line, Jf„ rcf. lLucvldl999h - a hat ( ) here denotes an intensity 
in the co-moving frame, 6 is for "blue" and lu indicates the tran- 
sition between some lower level I and upper level u\. The second 
one is the frequency-dependent, angle-averaged intensity j„. The 
procedure to obtain J v has been newly implemented; some details 
on this are given in Appendix [A] In order to reduce the noise on 
MC estimators, we simulate 500000 radiation packets in each MC 
run for the present study. 

2.1.2 Determination of the gas state; integration of the NLTE 
module 

The MC estimators are used to obtain an improved gues s on the 
state o f the gas. In the codes o flMazzali& Lucvl d 19931) and lMazzalil 
d2000h . the excitation and ionisation state is calculated using a 
"modified nebular approximation". A radiation temperature Tr for 
each shell is computed such that the mean frequency of a black- 
body at Tr matches the value v estimated from the MC run, and 
the electron temperature T e is approximated as 

T e := 0.9 T R . (1) 

An equivalent dilution factor W is calculated from 

W ■ B(Tr) = J. (2) 

In the nebular approximation, W, Tr and T e completely determine 
the excitation and ionisation. 

Numerically, the solution for the excitation and ionisation 
state is obtained by iteration ("plasma state loop", Fig.Q}. The cal- 
culation in the nebular approximation consists of the steps (I), (Ha) 
and (V) in Fig. Q] The cycles are started with two initial guesses 
n<; [initialisation, replacing step (I)], the excitation/ionisation state 
is calculated, respectively, and a density of electrons n e ,free is ob- 
tained from that state. In general, for both ?i e guesses, n e ,f re e will 
differ from the guess. The "real" electron density of the plasma will 
always be in between an n e guess and the respective n ei f rcc . In or- 
der to numerically determine the "real" n e (for which n e = n ei f ree ) 
and thus the plasma state, we run the loop (I)-(IIa)-(V), where (I) 
is a Newton-Raphson (NR) step to make the function f(n e ) := 
n e — n e ,frcc vanish (with a final accuracy of 0.005 x n e ). The NLTE 
module has been integrated into this plasma state loop to calculate 
the ionisation/excitation state of selected species in NLTE. Every 
time when it is called, the NLTE module obtains the necessary data 
from the MC code [Fig. Q~| step (lib)] and external files and cal- 
culates the required occupation numbers / ion fractions [steps (III) 
and (IV) -see Sec. [272] . 

2.1.3 Global temperature convergence and output spectrum 

MC runs and updates of the plasma state for all cells are iterated 
in turn. This also involves updates of T p h to enforce the given SN 
luminosity Lboi by compensating for the actual amount of backscat- 
tering (reabsorption at the photosphere). After convergence of the 
temperatures Tr and T e has been achieved (i.e. the relative changes 
with each iteration are < 1% in every cell), the electron densities 

2 This discretisation only applies to the variables describing the state of the 
gas; the positions of photon packets and lines are treated continuously. 
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Figure 1. Plasma state loop of the spectrum synthesis code; integration of the NLTE module. The cycle is initialized by inserting two initial guesses for n e 
and executing steps (II) - (V) for both guesses. As soon as two values for n e ,f rec have been calculated, a Newton-Raphson (NR) method can be used to obtain 
better guesses for n e [step (I)], and the loop is repeated until convergence (see text). Note that there are two independent NR routines in the code - the "outer" 
NR [step (I)] and the "inner" NR solving the rate equations [step (IV)] . 



and the occupation numbers do not change any more. The output 
spectrum is then obt ained from a formal solution of the transfer 
equation dLucvll 199*91) . 



2.2 The NLTE module 

We now describe the NLTE module. First we discuss the necessary 
external input data (model atoms, non-thermal excitation rates), and 
then we explain the actual solver [Fig.[T] steps (III) and (IV)]. The 
implementation follows closely that of ILucvI dl99 ll) . but uses up- 
dated atomic models and numerical techniques. 



2.2.1 Model atoms for H and He 

Atomic data for H are relatively easy or obtain. We constructed 
a simple model consisting of 13 energy levels. Each level corre- 
sponds to a principal quantum number; substates are not taken into 
account separately. For all calculations, the atoms are assumed to 
be effectively redistributed among the substates according to the 
respective g-f actors by collisi ons. Rate coefficients for our 13-level 
atom are from lMaurerl fcOld) . except for collisional ionisation and 
excitation cross sections u sed to calculate non-th ermal and thermal 
rates, which we took from ljanev & Smith! ( fl993h . 

Our He I model consists of 29 levels with principal quantum 
numbers n ^ 5, plus eight higher levels. For n ^ 5, substates dif- 
fering only in the magnetic quantum number J are treated as one 
"effective" state (as described above). For each n with 6 ^ n ^ 8, 
we combine all singlet states into one effective state, and all triplet 
states into another one. We finally take into account an = 9 triplet 
and singlet state, respectively, which we assume to be in LTE with 
He II. Energy levels for He I have been imported from the NIST 
Atomic Spectra Database (see acknowledgements). 

We have compil ed radi ati ve excitation ra te s for He I 

from 



Drake & MortorJ d2007t) . lLach & Pachuckil fcOOlh and 
Thakkarl feOOa). Photoionisation cross sections have been 



ICann & 

taken from iHummer & Storevl l ll998l) . and in some cases from 
iFernlev et al" 1 1987b . For c ollisional excitation cross sections we 
use iRalchenko et all J2008h for all transitions between states with 



n ^ 4. When calculating non-thermal collision rates (cf. Sec. 
12.2.21 ) to states with n > 4, we estimate cross sections with 
their recommended scaling formulae. Rate coefficients for ther- 
mal collisional excitation to or between levels with n > 4, how- 
ever, have conveniently been tak en from the rate compilations 
CHIA NTI dDereetalj[l997l.l2009l) and TMAD dRauch & Deetier] 
120031) . Collisional ionisation cross sections have been obtained 
from lRalchenko et al .1 d2008l ) (n ^ 5; again with scaling formulae 
for n = 5); for thermally averaged rates from n > 5 we used the 
formulae of lMihalas & Stone] dl968h . All de-excitation and recom- 
bination coeffi cients are calcu lated from the respective upwards co- 
efficients (e.g. lMihalaslll978l) . 

He II is approximately treated as one state; He III is not con- 
sidered in our calculations. This is sufficient for all but the outer- 
most parts of our atmospheres; we demonstrated this explicitly in 
iHachingerlfcoilh . 



2.2.2 Deposition of the gamma-ray energy 

Besides the usual radiative and (thermal) collisional bound-bound 
and bound-free processes, atoms in SNe undergo excitation and 
ionisation by fast Compton electrons. These are produced with en- 
ergies well above 1 keV and up to a few MeV by gamma radiation. 
In order to take the effects into account, the kinetic energy gain of 
the local electron population per unit time and volume (gamma- 
ray deposition rate or Compton heating rate He) must be known. 
Therefore, we model the bolometric light curv e simultaneously 
with the spectra. The light-curve code we use dCappellaro et al .1 
Il997lfl simulates the creation and transport of the gamma rays and 
writes out the relevant data. 

The kinetic energy of the fast electrons is assumed to be de- 
posited locally, causing atomic excitation and ionisation processes. 



3 The code version used here h as had some updates, most importantly to 
the opacity description [similar to Mazzali et al. ( 2 00 ll) . wit h a moderate de- 
penden cy on a temperature estimate, following the ideas of lKhokhlov et alj 
1 19931)1. 
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In order to calculate the resulting upwards rates, we have to deter- 
mine what fraction (deposition fraction) of the heating power 
He is deposited in each channel k, corresponding to an excita- 
tion/ionisation process or to Coulomb scattering with thermal elec- 
trons. This is achieved by simulating the downscattering of the fast 
electrons (see Appendix IB! for different plasma compositions. The 
results have been tabulated and are read in by the NLTE module, 
which calculates the energy H k (per unit time and volume) spent 
for process k to happen: 

H k =H c x D+, 
and uses this to compute the respective rates (see next Section). 



2.2.3 Solution of the effective rate equations 



In order to obtain the excitation/ionisation state of the elements 
treated in NLTE, we solve the rate equations in statistical equilib- 
rium. We abbreviate them here as: 



-4n = y 



(3) 



where A is the rate matrix and the n contains the number densities 
of atoms/ions state by state.The rate coefficients in A have (ther- 
mal) collisional and radiative contributions. The vector y contains 
the non-thermal collisional rates per unit time and volume. Only 
non-thermal rates from the ground state of an atom/ion are consid- 
ered significant (cf. Appendix B). One (redundant) rate equation r 
per chemical species "S" is replaced by an equation of species con- 
servation, i.e. A r j = 1 for all states j in "S", and y r — ns with the 
usual notation. 

The coefficients of A (except for those replaced by species 
conservation coefficients) are finally calculated as: 



Aij = IZji + Cji for j / i, and 



(4) 



3& 



where dj are thermal collisional rate coefficients (as usual), and 
IZij are radiative rate coefficients (see below). The indices i and j 
denote atomic states (including ionised states). 

Some details in which our rate coefficients differ from the 

usual form shall be briefly mentioned. 

For radiative bound-bound processes we use [cf. e.g. iLucvl 



J2002I) , equations (20-22)] the effective rate coefficients [s 
atom] : 



Ti.lu.sf! 



Am 



m R fb 

-Din -Jlu 

9u 



P 



per 



(5) 



— B iu J lu p. 



In this notation, the upwards rate TZiu^ni only depends on the 
number density ni (in the lower state V) and the downwards 
rate only on n u , if the escape/penetration probability /3 (e.g. 
lLamers &Cassinellilll999l) can be assumed constant (e.g. in the 
limit of very low line strength - otherwise /3 depends on occupation 
numbers). 

Furthermore, we have implemented a special treatment of ra- 
diative recombination to the ground state, for which the continuum 
is often optically thic k so that reco mbinations are followed by im- 
mediate re-ionisation. ILucvl J 1 99 lh took this into account by setting 
the respective recombination rate to zero. We have attempted an 
approximate "next-order" treatment here, in which we multiply the 



recombination rate in each cell with an estimate of the probabilitjQ 
for the respective photons to leave the cell or to be absorbed in lines 
(i.e. for no re-ionisation to occur). Our method is quite rough, but it 
is sufficient: it turns out that in the relevant cells the recombi nation 
rate to the ground state is normally negligible, as assumed bv lLucvl 
Jl99lh . 

The elements of y (again except those used for species equa- 
tions) represent non-thermal excitation processes. For excited states 
k, the entries y k are calculated as: 

Vk = H k /W k , 

where W k is the work for the excitation process. Each y k is a rate 
from the respective ground state, which depends on the mix of 
atom/ion species and free electrons. The non-thermal rate from/to 
each ground state "1" of an atom/ion is determined by excitations 
and ionisations taking atoms/ions away, and by ionisations bring- 
ing less ionised atoms/ions to the ground state (except for the case 
of neutral atoms - here yi is always negative): 



2/1 



Hi, 



Hh 



Wei. £2 w * 



E 



fcmax is the number of states in the ion considered, and all variables 
refer to that ion - except for H mn i s and W[ m S which refer to the 

lower ion. 

Charge exchange reactions (e.g. ISwartj 1 1994) and con- 
tinuum destruction proce s ses for photon s emitted by lines (e.g. 
iHummer & Rvbic kl 1 19851 : IChugail)l987ah are presently not con- 
sidered in our calculations, and distinct line overlap phenomena 
are not treated as we use the common Sobolev approximation. We 
assume that the rates we consider are dominant and the processes 
neglected have only a small impact on the results. Wherever rele- 
vant, we assume thermal electrons to have T e = 0.9 Tr as in the 
MC code; tests have shown that our results are not very sensitive to 
this choice. 

In Eq. [3] both y and A (via the escape probabilities) depend 
on the number densities n. Therefore the rate equatio ns are strong ly 
nonlinear, and we solve them using a NR method fas lLucvl dl99lh l. 
First, we test thermal ionisation/excitation patterns at different tem- 
peratures, and take as an initial guess riant the one for which 
\\A n — y|| (euclidean norm) is minimum [Fig.Q] step (III)]. Then 
we solve the rate equations [step (IV)] using the NR s olver package 
NLEQ2 dNowak & Weimanrj|l99ll ; lDeuflhardll26ol) . This imple- 
mentation includes sophisticated adaptive damping, and thus has an 
enlarged convergence radius; furthermore the solver is optimised to 
handle ill-conditioned Jacobians. When the accuracy estimator of 
NLEQ2 indicates convergence to 0.01%, we accept the solution. 



3 MODELS FOR SNE 2008AX AND 19941 

We begin by presenting models fo r well-observed reference ob- 
jects: the Type lib SN 2008ax (e.g. iTaubenberger et alj 201 1 ) an d 
the narrow-lined Type Ic SN 19941 (e.g. iFilippenko et al.lll995[) . 



4 To obtain this estimate, we calculate optical depth values in the contin- 
uum along rays, and finally average the probabilities p = e~ T over differ- 
ent directions and emission frequencies. For each ray considered, the optical 
depth integration ends (i.e. the photons are considered to have "escaped") 
either when the edge of the computational cell in the MC code is reached, 
or when a large line optical depth (~ 5) has been encountered, or when 
the photons have Doppler-redshifted out of the continuum in the co-moving 
frame - whichever happens first. 
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The observations we use have been carried out at epochs between 
16.0 d and 40.6 d past explosion. Thus, they cover the epochs in 
which the H and He lines in SNe Ilb/Ib are usually strongest and 
can be used best for abundance determination. 

The initial ejecta models we assume for the two SNe have 
abundances and densities resulting from 1-D hydrodynamics/nu- 
cleosynthesis calculations. Specifically, for SN 2 008ax we use 
the explosion model 4H47 ( Shigevam a et al.|[l994l) constructed to 
explain a relatively similar SN lib (SN 1993 J). For SN 19941, 
llwamoto et al ] dl994h have computed the explosion model C021. 
These models may not be the most recent ones, but they have a 
structure consistent with the observations. 4H47 and C021 actu- 
ally use the same stellar core; only the mass and the composition 
of the envelope is different. The kinetic energy of the ejecta in the 
models considered here is always ~ 10 51 erg. The 4H47 model we 
use has been somewhat smoothed so as to avoid sharp line compo- 
nents due to overdense narrow shells in the 1-D model. 

Starting from the initial values of the nucleosynthesis models, 
we fit the time series of spectra u s ing th e abundance tomography 
technique devised by IStehle et alj d2005l) . The method constrains 
the abundances in the outer layers of the SN starting from the ear- 
liest observations, and mapping the composition of deeper layers 
from spectra observed later (when more of the ejecta is optically 
thin). 

In this work, we search for best-fit abundances under the con- 
straint that the stratification typical of a core-collapse SN shall be 
preserved. This means that increasingly heavier elements dominate 
the ejecta from the outside to the inside, and that only species abun- 
dant in typical core-collapse ejecta (H, He, C, O, Si, Fe-group) are 
allowed to be dominant. A basic tomography is conducted with one 
abundance zone per observed spectrum. In the present models, we 
augmented the number of zones by adding two zones on the out- 
side (i.e. the earliest spectrum was modelled with three abundance 
zones). This significantly improves so me of the spectral fits (see 
below). We ran the light curve code o f lCappellaro et all d 19971) on 
the current ejecta model from time to time as the spectral modelling 
process went on, in order to ensure compatibility with the observed 
light curve and to make sure the non-thermal excitation rates are 
consistent with the ejecta model. 

The density profile of the respective explosion model was left 
unchanged in the modelling process. Probing the large parame- 
ter space opening up otherwise would have been computationally 
costly and is not in the focus of the current work. 

Below, we show and discuss the models for the two SNe. The 
code input parameters of the models are shown in Appendix Icl and 
the density profiles are shown in Section|4] 

3. 1 SN 2008ax / SN lib model 

SN 2008ax has been extensively observed by several groups 



dPastorello et al 



been extensively observ ed by several g roups 
20081; iRoming et alj|2009l IChornock et alfeoict 



iTaubenberger et alj201 1 ). The observ ations used for ou r modelling 



are combined from all these sources faachingenl20l 1) to obtain a 
good wavelength coverage at four epochs [16.0 d, 22.1 d, 30.6 d 
and 40.6 d after t he explosion, which we assume to have happened 
at JD 2454528.8 dPastorello et alj|2008l) l. In order to homogenize 
the data, we multiplied spectra which did not match contemporane- 
ous photometry with simple correction functions /(A) (with values 
close to unity). For all synthetic photometry, we used IRAF with 
SYNPHOT and TABLES (see acknowledgements). 

We assume a total reddening of E(B — V^tot = 0.30 mag 
(with a Galactic contribution of 0.02 mag) for the SN, and a dis- 



tance modulus of (i = 29.92 mag [all values from lPastorello etHll 
<f2008h l. 

The spectral models for SN 2008ax are shown in Figure|2]with 
the observed spectra and line identifications. All data in the plot are 
de-reddened and de-redshifted. 



3.1.1 t — 16.0 d spectrum 

At this epoch, the dominant features are made by He and H, apart 
from some Ca, Mg, Si and Fe lines in the B and V bands. The He I 
lines are relatively weak, as significant amounts of He are below 
the photosphere (v p h = 8400 km s _1 ). Also, the radiation temper- 
ature is generally very high, which leads to high ionisation rates 
from excited states of He I, decreasing the occupation numbers of 
these states. Yet, numerous He I lines can be identified (Fig. [2] red 
marks). The features at A6678 A7065 and AA20581, 21120, 21132 
allow for an unambiguous detection of He, especially when ob- 
served together. Other features may be contaminated by lines of 
different elements more frequently (see Sec. 13.1.31 . 

We inserted two additional abundance zones above the pho- 
tospheric one, with lower boundaries at 10800 and 13000 km s _1 , 
respectively. This allows for a better match to specific spectral fea- 
tures. In the zone between 10800 and 13000 km s" 1 , the abun- 
dances of C, O, intermediate-mass elements (IME: Ne...Ca) and 
Ni decrease strongly; this is needed in order to avoid high- 
velocity IME lines (especially Ca) and excessive line blocking (as 
caused by 56 Ni). In the outermost zone, these abundances are fur- 
ther reduced and replaced by H, which makes the high-velocity 
wings of the H lines a bit stronger. 

The model matches the observations well in the optical and 
the IR, except for the fact that H lines are generally somewhat 
too strong and He lines a little too weak. We accept such small 
mismatches here, as they depend on details of the model (e.g. the 
gamma-ray flux due to 56 Ni or the density distribution). Around 
3500 A the flux of the model is somewhat too low which indicates 
too much absorption by Fe-group elements. The abundances we as- 
sumed in the model were however necessary to fit the light curve 
reasonably well. 



3.1.2 t = 22.1 d spectrum 

At 22. 1 d after explosion, the SN is 3 d past its B maximum. The 
spectrum is cooler than before and spectral features get deeper, as 
more material is uncovered (v P h = 7300 km s _1 ). 

The high-velocity wing of the strong He I A10830 feature, 
which is too deep in the model, is produced in the outermost regions 
of the ejecta. In these regions, the density in the model is proba- 
bly too high. With a lower density, the material in the outer lay- 
ers would not only be more dilute, but also may become extremely 
ionised (up to He III, which we do not consider here) at the expense 
of He I, The early-time light-curve spike of SN 2008ax is much 
narrower than that of SN 1993J. This also i ndicates that SN 20 08ax 
had an exceptionally low envelope mass (cf. lArcavi et al .1201 ll) . i.e. 
lower densities than in the 4H47 model, which originally has been 
constructed to explain SN 1993J (Shigevama et al.lll994). 



3.1.3 t — 30.6 d spectrum 

By day 30.6, the photosphere has receded to v p h = 6000 km s _1 . 
In general the spectrum shows deeper lines than those at earlier 
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Figure 2. Spectral models for SN 2008ax (red lines) compared to the observations at four epochs (black lines). The insets show the infrared (IR). Line 
identifications are given for prominent lines at the first and third epoch - red: He I lines, blue: H lines. The H« line first remains constant and finally gets a bit 
weaker with time, while most He lines gain additional strength. 
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Figure 3. Bolometric light curve calculated from our ejecta model of 
SN 2 0O8ax compared to ob s ervations [based on compiled UV-to-IR data 
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Figure 4. Abundance tomography plot for SN 2008ax, as derived from the 
fit to the spectra and the light curve. 



epochs, and a multitude of strong He I features is visible. The red- 
der colour and lower flux level of the spectrum lead to a lower ion- 
isation. This makes the occupation numbers in the excited states of 
He I higher, with the effect of stronger lines. The clearest He fea- 
tures are at AA6678, 7065, 7281 and around 20000 A. The A 10830 
featur e may possibly have a contribution from C I (Sa uer et"al] 
l2006h ; however in our simulation the lower levels of the respec- 
tive lines (C I AA10691, 10730) are not sufficiently populated. The 
A5876 feature is extremely prominent in the optical, but has some 
contribution of N a I D, which is common in core-collapse SNe (cf. 
ISauer et alj|2006h . 



3.1.4 t = 40.6d spectrum 

40.6 days past explosion, the development described above con- 
tinues and the He lines reach their maximum strength (compare 
the He I A6678 line with Ha in Fig. [2] lower panels). The photo- 
sphere (v p h — 3100 km s~ ) has reached the bottom of the He- 
dominated zone in the hydrodynamics/nucleosynthesis models of 
dShigevama et aljl994l) . Despite the late epoch, the observed spec- 
trum is quite well matched. The match to the IR spectrum is par- 
ticularly good, although the models were crafted so as to primarily 
optimise the fit to the optical observations. 

The IR spectrum now shows a lot of structure, not only with 
the marked He I feature at ~ 20000 A, but also with other features 
due to He, Co, Mg and more elements. 



3.1.5 Light curve 

Our model light curve is shown together with the bolometric LC of 
SN 2008ax in Fig. [3] The match is generally very good, taking into 
account that the opacity description used for the calculation is only 
approximate. The total 56 Ni mass in our model is O.O9M0. 

The earliest points of the observed curve might be better re- 
produced if the model envelope was somewhat less massive (cf. 
Sec. 13.1.21 , so that radiation could escape faster. In addition, these 
points may have a contribution from the envelope heated during the 
explosion (as in Type II-P/L SNe), which we neglect in the model. 



3.1.6 Abundance stratification 

The abundances in our ejecta model are plotted in Fig. [4] With a 
photospheric velocity of 3100 km s^ 1 at +40.6 d, our tomography 
only samples the He-dominated part of the SN. The total mass of 
He in the model is ~ 1.2Mq. 

In the outer envelope (v > 7500 km s _1 ), the H mass fraction 
is ~ 10-40%. A larger mass fraction of H would lead to the forma- 
tion of a deeper Ha line, which is already somewhat too strong in 
the models. The total H mass in our model is only 0.05 A/ Q , which 
is less than in the original 4H47 ejecta model. However, this is in 
line with another analysis of SN 2008ax dTsvetkov et al.l2009h . The 
remaining mass fraction has been attributed to He, which at these 
velocities does not contribute to line formation (with the A 10830 
high- velocity feature as an exception). The actual H/He envelope 
of SN 2008ax may indeed be somewhat less massive than in our 
model and have a larger H mass fraction instead. 

The fraction of Fe-group elements in the outermost part of the 
ejecta is sub-solar (~ | of the solar values). This reduces the flux 
blocking in the outer layers with respect to the solar-metallicity 
case. In a solar-metallicity model, line blocking would cause more 
backwarming (a more intense, "hot" radiation field within the at- 
mosphere), which leads to a somewhat altered ionisation state and 
a worse spectral fit. Not only H and Fe-group elements have been 
reduced in the outer envelope: for Ca, a reduction to less than 10% 
of the solar value was required in order not to produce spurious 
high-velocity absorption. All this again indicates that the density 
(and mass) of the outermost layers in the 4H47 model is higher 
than in the observed object. 

Most line strengths are determined by the abundances in the 
intermediate part of the ejecta (v ~ 4000- 12000 km s _1 ). Many 
elements have an influence on the spectra in this velocity range, 
except for Si and Ne whose abundances are therefore estimates. 

From the model, it is evident that the composition is relatively 
strongly mixed, with some percent of IME and Fe-group elements 
almost everywhere. Some outward-mixing of 56 Ni has to occur for 
the light curve to rise fast enough; we have assumed this to be as 
strong as allowed by the spectra. 

3.2 SN 19941 / SN Ic model 

SN 19941 is a v e ry well-obs e rved l ow-mass Type I c SN 
dYokoo et al.1 1 1994 | Sasak i et al ] 1 19941 : iFilippenko etai] 1 19951 : 
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Figure 5. Spectral models for SN 19941 (red lines). Observations or spectra interpolated from observations at adjacent epochs (cf. text) are shown in black. In 
the insets, the IR is plotted (no observations available). Line identifications are given for prominent lines at the first and third epoch - red: He I A10833 line. 
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IClocchiatti et alJll996l ; lRichmond et"aTlll996h . Thanks to the good 
time coverage and quality of the observatio ns it has already been 
subject to severa l radiative transfer studies dlwamoto et al.|[l994l ; 
ISauer et al .120061 and references t herein ). We use the spectroscopic 
observatio ns bylff ilippenko et al .Id 19950 and the photometric obser- 
vations bv lRichmond et al .1 119961) , We calculate model spectra for 
SN 19941 at four epochs. In order to render the models in this paper 
well comparable to one another, the epochs were chosen to match 
exactly those of SN 2008ax (we a s sume that SN 19941 exploded on 
JD2449438.fJJ)- iFilippenko et alj l l 19951) have taken a spectrum at 
an epoch 16 d, which we fit with our first model spectrum. For the 
subsequent epo chs, we compare o ur models to spectra interpolated 
[as described in lHachingeJ HoTlh l from observations taken sooner 
and later. The first spectrum and the three interpolated spectra are 
again calibrated against contemporaneous photometry. IR spectra 
are not availabl e for SN 19941. 

Following TSauer et al.1 d2006h . we assume for SN 19941 a to- 
tal reddening of E(B — V) to t = 0.3 mag [the Galactic reddening 
is 0.035 mag, ISchlegel et ail dl998h l. and a distance modulus of 
H = 29.60 mag. However, all models and observed spectra plotted 
below are rescaled to the distance of SN 2008ax (p — 29.92 mag) 
and printed on the same scale in order to facilitate visual compar- 
ison. Again, all spectra have been de-reddened and de-redshifted 
prior to plotting. 

The spectral models for SN 19941 with line identifications are 
shown in Fig. [5] 

3.2.1 t = 16.0 d spectrum 

The SN spectrum at t — 16.0 d, which is already ~ 5 d past max- 
imum for this SN and corresponds to a photospheric velocity of 
8900 km s~ , shows a mix of O, IME and Fe-group lines. In the 
region around 4300 A, the SN shows a deep trough (which is due 
to Mg II, Cr II and Ti II) similar to 91bg-like SNe la. Some hint of 
a C II line is visible at the red edge of the weak Si II A6355 line. 

Model line velocities in general tend to be a bit too high, which 
reflects the fact that the decline of the density of model C021 in the 
outer layers is extremely shallow. With respect to the Si II A 6355 
feature, where the mismatch is largest, iBranch et al. I d2006l) have 
furthermore suggested a significant contribution of H mixed into 
the star's core. We have assumed H not to be present in our models. 

The model includes a tiny amount of He in the intermediate 
and outer layers (0.04M Q in total); this makes the He I A10830 
line appear, but none of the weaker He I lines. The observed spectra 
at later phases indicate the possible presence of the A10830 line; 
unfortunate l y, the feature is never fully covered by the observations. 
Saue r et al. I d2006h showed that the observations in this wavelength 
range may also be reproduced by C I lines instead of He I, if the 
excitation/ionisation state of C is appropriate. 

Again, we have used two abundance shells above the 
t — 16.0 d photosphere (with lower boundaries at 15400 and 18500 
km s , respectively) to optimise the fit to the spectrum. In these 
zones we have assumed that the composition is dominated by C 
and contains some He, as was the cas e for the outer shells in the 
C021 nucleosynthesis calculation of llwamoto et alj l ll994l) . The 
density of IME and Fe-group elements is reduced as in the model 
for SN 2008ax. This has merely been done to maintain a consistent 
approach (e.g. the Fe-group elements are also reduced down to ~ | 

5 This date is obtained subtracting a rise time of 12 d jl wamoto et al 1 19941) 
from the date of B maximum, JD2449450.56 iRichmond et all 1996ft . 
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Figure 6. Bolometric light curve calculated fro m our ejecta mod el of 
SN 19941 compared to observations (as shown by ISauer et alj|2006l with 
his IR corrections based on SN 2002ap). The 56 Ni mass is O.O7M . 

of the solar abundance); the exact amount of the reduction plays 
less of a role for the SN 19941 models. Only the Ca mass fraction 
has to be very much reduced in the outer layers (to X (Ca) ~ 10 -6 ) 
in order not to produce high-velocity features. 

3.2.2 t = 22.1 d spectrum 

The 22.1 d spectrum of SN 19941 is significantly cooler (redder), 
as the luminosity of the SN has declined rapidly. The photospheric 
velocity has decreased to the low value of v P h — 3100 km s , 
which is plausible considering the low ejecta mass of SN 19941. 

The IR spectrum begins to show a lot of structure owing to 
lines of Co II and Mg II. 

3.2.3 t = 30.6 d and t = 40.6 d spectra 

At 30.6 d, the photosphere has receded to 1400 km s~ . The limit 
of applicability of our code is certainly reached here, as a major 
fraction of 56 Ni is now above the photosphere, but nonetheless we 
have been able to calculate a decent model spectrum. We continued 
the calculations to t = 40.6 d (v p h = 910 km s _1 ). Other models 
in the sequence from SN 2008ax to SN 19941, which we construct 
below, do not have such a deep-lying photosphere at these epochs 
and therefore will be more reliable. 

The composition near the photosphere at these epochs is com- 
pletely dominated by 56 Ni, which is required to make the light 
curve sufficiently bright. The sensitivity of the optical spectrum on 
the composition in the photospheric layers is limited. Instead, prac- 
tically all elements abundant in the intermediate and also in the 
outer region of the ejecta contribute to line formation. 

The He in the outer layers of the model leaves a very broad 
A 10830 line now, with a high-velocity component a bit too pro- 
nounced. 

3.2.4 Light curve 

The light curve calculated from our SN 19941 model (Fig. [6jl pro- 
vides a decent match to the observations in the early phase. After- 
wards, at 30 d past maximum and later, there is a mismatch between 
model and observations. We could not increase the 56 Ni mass in the 
model (0.07A/q) to improve the match because the light curve peak 
is already a bit too high. 

The mismatch between peak and tail indicates that the model 
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Figure 7. Abundance tomography plot for our SN 19941 model. 



is not opaque enough, although we have already concentrated the 
56 Ni in the centre of the ejecta. The model does not delay the es- 
cape of photons to late enough epochs. However, in models with 
moderate differences to ours in the opacit y description and/o r the 
ejecta configuration dlwamoto et alj Il994l : ISauer et al.l l2006h . the 
light curve tail has been successfully matched. Thus, we can as- 
sume that the density and chemical composition of our model is all 
in all reasonable. 



3.2.5 Abundance stratification 

The abundances that we obtained from fitting SN 19941 (Fig. [7) 
can be regarded as typical for SNe Ic. Helium appears as a trace 
element in the outer layers (total He mass: 0.04 Mq); these layers 
are dominated by C so that C II lines can form. A hint of the C II 
lines, which are rather weak intrinsically (i.e. in terms of oscillator 
strength), is visible redwards of Si II A6355 in the spectra. Below, 
O is the dominant element; mass fractions of IME and Fe-group 
elements increase further inwards. 

The light curve requires a strong concentration of 56 Ni and 
other Fe-group elements in the core. Accor ding to our calculat ions, 
this is compatible with the spectra, while ISauer et al] d2006l) had 
to reduce the Fe and 56 Ni mass fractions in their late-time models. 
This reflects the fact that lSauer et alj d2006h used one-zone models 
(i.e. homogeneous abundances throughout the envelope), and that 
too high a mass of Fe-group elements in the ejecta is not allowed 
by the spectra - even if the material becomes more dilute at late 
epochs. 
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Figure 8. Energy rescaling function r)(v) used to morph the "core" 
(Fe,Si,0,C and He zones) of the 4H47 model into the shape of the M4.0 
model (upper panel). Lower panel: same for the M4.0 — > C021 transfor- 
mation (where the core consists of the Fe.Si.O and C layers). 

4.1 Set-up 

The model sequence we set up is itself based on progeni- 
tor/explosion models which have the same stellar Fe/Si/O/C core 
and only differ in the envelope. We have decided to present iso- 
energetic models (i.e. the explosion energy is always 10 51 erg), 
which matches our choice of reference objects (SNe 2008ax and 
19941). The density/abundance structure of models in the sequence 
is set up so as to make a smooth transition between the models for 
SNe 2008ax and 19941. 



4 MODEL SEQUENCE - FROM SNE IIB TO SNE IC 

We now present our model sequence representing SNe with differ- 
ent degrees of envelope stripping. First, we explain how we obtain 
density profiles and the code input parameters for the models. Then, 
we show the most interesting synthetic spectra from the sequence 
and estimate the observable consequences of the different He (and 
H) content. The sequence has SN 2008ax and SN 19941 as start 
and end points, respectively. It is helpful for the construction of the 
sequence that the explosion models 4H47 and C021 used for these 
SNe are based on the same stellar core. 



4.1.1 Density profiles 

In the context of models with constant kinetic energy, the step from 
a SN with a massive envelope to a more stripped SN is mostly re- 
flected in the structure of the remaining H/He envelope and of the 
outer core. These will carry much more kinetic energy per unit mass 
in a more stripped SN. The innermost layers will show less of a dif- 
ference [in_their p re-explosion state a s well as their kinematics, cf. 
lArnettl dl974h and lMatzner & McKeel dl999h l. To construct our se- 
quence, we all in all have to perform a simultaneous reduction of 
the envelope mass and an adequate redistribution of the remain- 
ing ejecta in velocity. The first aim is to gradually transform the 
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Figure 9. Density profiles for ou r SN Ilb/Ib/Ic models IShigevama et alj 
ll990Ul994l ; llwamoto et a?1ll994) and two transitional cases. The transi- 
tional profiles are constructed as described in the text. Major abundance 
zones are indicated. Note that most of the exploded C/O core of the star 
will consist of O in all models according to the original nucleosynthesis 
calculations (and also according to our spectral models). 



4H47 (SN lib) ejecta model into th e "M4.0" (SN lb ) mode l pro- 
duced from the same stellar core bv lShigevamaetalJjl99Ch . Sec- 
ond, we transform the M4.0 density profile so that we obtain the 
SN Ic model C021 in the end. 

For each of these transformations, we first rescale the core 
of the ejecta (i.e. the ejecta without the envelope which is grad- 
ually to be removed). The relocation of the respective mass ele- 
ments in velocity space can be expressed in terms of an energy 
rescaling function 7](v) = 5r («) such that the kinetic energy of a 
mass element 8m at an original velocity v increases from E to E' 
during the rescaling process. We have (re-)constructed p iecewise 
linear rescaling functions 77 [Fig. [8j cf . lHachingeil d201 11)] for the 
two transformations "4H47 core" — > M4.0 and "M4.0 core" — ► 
C021. For each of the two transformations, these functions are ap- 
plied step by step, so as to obtain 20 models, respectively [i.e. an 
application of the function 77(77) = 1 + (r)(v) -l)x = brings us 
n steps out of 20 towards the next model]. After the core rescal- 
ing, we attach a newly-constructed, less massive envelope to each 
model. Depending on the respective core configuration, this enve- 
lope is chosen to make the final kinetic energy of every (complete) 



model equal to 10 erg. For the construction of the outermost part 
(high-veloc ity ejecta, v > 2000 km s -1 ) we use Equation (8) and 
Table (3) of iMatzner & McKed dl999r) . This part is then connected 
to the rescaled core using a simple density function of the form 
p(r) = a x exp(— r/b) + c, yielding a reasonably smooth tran- 
sition. The original density profiles and an example model within 
each sequence are shown in Fig. [9] The displacement of different 
zones as velocity is rescaled is also indicated in the Figure. Each se- 
quence model has been smoothed somewhat in order to avoid jags 
which result from small scale variations in the initial (4H47 / M4.0) 
density profiles and from the not entirely smooth 77(77) functions. 



4.1.2 Preparation of the SNe 2008ax and 19941 models and final 
set-up of the sequence 

In order to fully specify the model sequence, we now must define 
a sequence of abundance stratifications, and values for the photo- 
spheric velocity and luminosity at each epoch. 

To do this, we first identify "corresponding shells" in our 
ejecta models for SNe 2008ax and 19941. The ejecta of the two 
models have originally been divided into six shells with different 
abundances. By setting up both ejecta models to have eleven abun- 
dance zones, we have been able to establish a one-to-one correspon- 
dence, in the sense that the 71-th shell of the SNe 2008ax model 
starts and ends at the same enclosed mass coordinate as the n-th 
shell of the SN 19941 model (for shells at M < 0.9M©). In Ap- 
pendixlCl(Table lCTt , we give the parameters of the SN 2008ax and 
SN 19941 models zone by zone with some technical annotations. 

The abundance zones of the models within the transition se- 
quence are then constructed by moving the abundance zone bound- 
aries with the mass elements and interpolating the abundances in 
each zone between the SN 2008ax and SN 19941 values (linearly, 
in 20 + 20 = 40 steps for the entire transition). 

Finally, we have to define the luminosity and the photospheric 
velocity of the sequence models for each epoch. The luminosity 
is determined by linear interpolation, just as the abundances. The 
photospheric velocities are calculated, epoch by epoch, using the 
Stefan-Boltzmann law, inserting the new L value and a linearly- 
interpolated effective temperature T p h (again in between the val- 
ues for SN 2008ax and SN 19941). This procedure has the advan- 
tage that photospheric black-body temperatures change in a well- 
defined manner. Our photospheres thus defined for the sequence 
models do not necessarily correspond to the abundance shells con- 
structed before. Therefore, we just insert them as lower boundaries 
for the spectrum calculations, without any impact on the abundance 
structure. 



4.2 The SN lib 



SN lb transition 



Figure [Tol shows a selection of models representing the sequence 
from the SN lib model over two transitional models|j ("SN Ilb/Ib" 
and "SN Ib/IIb") to the SN lb model. The Ha line becomes pro- 
gressively weaker and the H/3 line, which is visible at the earliest 
epoch, vanishes (see marks in the Figure). 

Finally, the SN lb model still has residual absorption in the 
Ha region. This is due to Si II A6355 and a (weaker) "real" Ha 



B These models have been identified by eye as those for which a clear H 
identification is still possible ("SN Ilb/Ib", model no. 6 in the transition 
sequence — > 20) or already difficult ("SN Ilb/Ib", model no. 10). 
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Figure 10. SN lib — > SN lb sequence of spectral models. The model designations refer to the SN type the respective model represents (two types in case of 
transitional models). Line identifications are given for a few He and H lines with which the SN class can be well determined. In the uppermost panel, an inset 
shows the H/3 region in more detail, and we have also marked the Si II A6355 line which appears in the SN lb model, while it vanishes in the high-velocity 
wing of H« for more H-rich models. 
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Figure 11. SN lb — ► SN Ic sequence of spectral models. The model designations again refer to the represented SN type (two types for transitional models). 
Line identifications are again given for the He features most useful for classification. 
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absorption, as some mixing of H into the He-rich zone has been 
assumed. 

The H mass for which the Ha/Si feature becomes dominated 
by H and for which H/3 absorption appears is between 

M(H)sNiwnb = O.O25M 

and 

M(H) S Niib/ib = 0.033M Q 

in our model sequence. 

The values should be taken to refer to the outer envelope; a 
small mass fraction of H throughout the SN may not have any visi- 
ble effect. Furthermore, charge exchange may decrease the H I/H II 
fraction if other elements such as He are present, and differences 
in the 56 Ni distribution may change the H excitation and ionisa- 
tion. The point here is, however, that a relatively H-rich envelope 
[X(H) > 10%] with a very small H mass is sufficient to make the H 
lines sh ow up. This is consi stent at least with some progenitor mod- 
els (e.g. lGeorgy et al.l2009l) . and with observations: it has been sug- 
gested that a weak Ha lin e is actually present at least in some SN e 
classified as Type lb (e.g. lHamuv et alj|2002l. IValenti et alj|201ll) . 
Apart from the presence of Ha, there is not much difference be- 
tween SN lib and SN lb spectra some weeks past explosion. 

4.3 The SN lb ► SN Ic transition 

As the He is removed from the models (Fig. II lb . we see a very 
pronounced change in the spectra. C, O, intermediate-mass and Fe- 
group elements begin to exclusively shape the spectrum (cf. Sec. 
13.21 ). Depending on the actual composition and the individual shape 
of the envelope, the sequence models develop more or less pro- 
nounced high- velocity components in the Ca IR triplet and the He I 
A 10830 line. Although quite prominent, these are not really signif- 
icant for our study. We focus on the visibility of clear He features. 
Again, we consider two transitional models ("SN Ib/Ic" and "SN 
Ic/Ib" models in Fig. [TT] — these are models no. 23 and 28 in the 
sequence 20 — > 40, in which the He abundance declines rapidly at 
the beginning). 

While many He I lines are possibly blended with lines from 
other species, we have marked three features in Fig. [TT] (A6678, 
A7065 and AA20581, 21120, 21132) which are good indicators of 
He content, especially when taken together. The IR feature (A20581 
with minor contributions from AA21 120, 21 132) is particularly use- 
ful, as it is strong, but does not saturate too quickly. Before max- 
imum light, and for some weeks after maxi mum, all He I lines 
tend t o become more pronounced with time (cf. lTaubenberger et al.l 

I20T1I) . 

The He mass for which we begin to see the three features as 
pointed out above is between 

M(He)sNMb = O.O6M 

and 

M(He)sNiMc = 0.14M S . 



5 DISCUSSION 

Our model sequences have allowed us to obtain mass limits for H 
and He in low-mass core-collapse SNe. When He lines are not ob- 
served, the He content must be in the order of 0. 1 Mq or less. This is 



particularly interesting: stellar models have difficulties in explain- 
ing a strong loss of He, while a sma ll H mass can be reproduced 
dGeorgy et all2009tlYoon et alfcoioh . 

Unless 3-D effects play a very strong role, our He mass limit 
is quite robust (while H limits are more uncertain, cf. Sec. \A.2\ . 
We would like to further discuss this aspect and the state of He 
in typical models. Finally, we assess whether current progenitor 
models match our results, and discuss the importance of He I IR 
lines as fine diagnostics on the He content of the progenitor stars. 

5.1 State of He in our ejecta models; robustness of our results 

In order to test the reliability of our conclusions, we have calcu- 
lated the "SN Ib/Ic" model again with heating rates He artificially 
reduced by a factor of five throughout the model atmosphere. Re- 
markably enough, we found that the impact of this extreme change 
on our synthetic spectra is practically negligible. Any equivalent 
effect on the excitation/ionisation state of the He-rich layers will, 
therefore, only have a mild impact on our conclusions. 

To illustrate the reason for this, we have plotted occupation 
numbers for He in the "SN Ib/Ic" atmosphere at 22.1 d past ex- 
plosion (Fig. |T2] upper panel). This is an epoch at which the He 
lines are already strong, and at which observations are available for 
many SNe Ib/c. 

First, it is surprising that most of the He is actually singly 
ionised (Fig.[T2] upper panel, red curve). Also, the occupation num- 
bers in the excited states are large, with departure coefficients [i.e. 
occupation numbers with regard to Saha-Boltzmann numbers at the 
electron temperature T e (r)] in the order of 10 3 in deep atmospheric 
layers and 10 7 in the outermost layers. At 40 d past explosion, the 
departure coefficients have increased to 10 10 in most of the He-rich 
zone. 

The strong departures from LTE are due to the structure of the 
He atom and the typical physical conditions inside a SN. The be- 
haviour of He can best be understood as follows: The He I excited 
and ionised states are separated from the ground level by large en- 
ergies (e.g. lBashkin& StoneJI 19751) . In addition, transitions from 
the two lowest excited states of He I (2 3 S, 2 1 S) to the ground state 
are forbidden; transitions from the ionised state to the He I ground 
state are ineffective because of direct re-absorption. Therefore, the 
He I ground state is only sparsely coupled to all other states. On the 
other hand, the He II ions together with excited He I atoms form a 
strongly coupled system in which radiative transitions establish a 
rate equilibrium. The result is then determined by: 

(i) the equilibrium between He II and the excited states of He I. 
The energy differences within this system correspond to 
UVOIR wavelengths, at which the SN has a strong radia- 
tion field. The ratios of the respective occupation numbers 
are therefo re well estimated by th e nebular approximation 
(Sec. l2"X2irMazzali & Lucvll993l) . Owing to the fast radia- 
tive rates, this equilibrium is essentially undisturbed by the 
transitions to and from the He I ground state in the present 
model. Only in few cases (depending on model and shell), 
the downwards rate from some of the excited states can be 
large enough for the occupation number in these states to 
deviate significantly from the nebular one. 

(ii) the total upwards rate from the He I ground state to the ex- 
cited/ionised system. This rate is practically only due to fast 
electrons, as the energy gaps to the excited states are very 
large, and radiative transitions from the ground state to the 
lowermost excited states are forbidden. 
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Figure 12. State of He in a typical model, as a function of velocity (radius). 
Upper panel: SN Ib/Ic model at 22. 1 d - normal version; lower panel: com- 
parison to the modified model version with Compton heating reduced by a 
factor of five (thick, lighter lines). We have plotted the number densities of 
He I atoms in different states (ground state "GS"; states 2 3 S, 2 X S, 2 3 P, and 
2 P), as well as the number densities of He II in the GS, of all atoms/ions 
in total (n to i - regardless of NLTE/nebular treatment), and of free electrons 
(n e ). The step in the quantities at ~ 6000 km s — 1 corresponds to a change 
in the He number fraction from 22% to 1 % in the model. 



(iii) the total downwards rate from all excited/ionised states 
to the He I ground state. The occupation ratio of the ex- 
cited/ionised system to the He I ground state is determined 
by the equilibrium of the upwards rate (ii) with the down- 
wards rate. 

In the lower panel of Fig. [12] we show what happens if the 
original model (upper panel) is modified reducing the Compton 
heating rate Ho as discussed. The ratio of He I atoms to excited 
and ionised atoms changes. However, as long as He is mostly singly 
ionised, this only has a small influence on the number of He I atoms 
in excited states, which is in equilibrium with He II via radiative 
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transitions. As long as He II is sufficiently dominant, the optical 
depths of the He lines in our spectra, which all originate from ex- 
cited He I states, will be less dependent on the exact heating rate, 
but mostly on the number density of He and the radiation field. 

The SN Ib/Ic model we used contains a moderate He mass 
fraction of 22% in most parts of the envelope. This corresponds to 
the assumption of relatively strong mixing, which is expected to 
make nonthermal excitation of He weaker as more energy of the 
fast electrons is deposited into other species. Our upper limits on 
the He mass are therefore relatively conservative. 



5.2 3-D effects and our results 

3-D effects in the ejecta configuration may play a signifi- 
cant role in core-collapse SN spectroscopy. Although extended 
H envelopes supposedly remain symmetric i n the explosion 
dLeonard & Filippenkol2005llLeonard et al.l2006t) . the He layers in 
the ejecta can have a thickness varying more or less with viewing 
angle. The heating due to radioactivity may be viewing-angle de- 
pendent as well. A larger He mass may be hidden if marked asym- 
metries are present. 

As we cannot precisely judge the possible asphericity from 
ab-initio explosion models yet, we just mention these facts as a 
caveat to our results. A quantitative assessment of viewing-angle 
dependencies based on a specific explosion model and 3-D radia- 
tive transfer treatment may be an interesting topic for a future study. 



5.3 Compatibility of our results with progenitor models 

Current stellar evolution models for c ore-collapse SN progenitors 
dGeorgv et al]|2009l ; pYoon et al.ll2010l) . at least those involving bi- 
narity, succeed in explaining envelope configurations similar to our 
SN lib and SN lb models. 

Reconciling our limit on He masses in SNe Ic (~0.06- 
O.14A/0) with the progen itor models is , howe ver, more difficult. 
The models presented by iGeorgy et al.l d2009l) . which are single 
Wolf-Rayet stars, never have less than 0.3Mq of He. In order to 
explain the total SN Ic rate with their progenitors, th ese authors 
need t o allow for He masses of 0.6Mq in SNe Ic. Also lYoon et al.l 
who simulated binaries with one to three phases of "case 
A" or "case B" mass transfer, produce only a small amount of pro- 
genitors with M(He) < 0.5 Mq. The models of both groups are, 
however, also too massive to explain SNe Ic as we studied them 
here. 



5.4 IR He 1 lines as diagnostics for the SN lb — ► SN Ic 
transition 

For extended studies on the H/He content of stripped-envelope 
core-collapse SNe, it is important to know which lines or spectral 
regions allow for an accurate measurement of abundances. For H, 
the Ha line is the best choice especially when abundances are low, 
but for He the question is more difficult to answer. 

There is a multitude of He I lines, but the most interesting fea- 
ture for diagnosing the He content is certainly the IR 2p one. The 
intrinsic absorption strength of the transitions, together with typical 
population of the lower levels, is obviously favourable. Therefore, 
small amounts of He (like those occurring in our SN lb — > SN Ic 
sequence) suffice to make the feature show up, but also for larger 
He masses the lines do not saturate too quickly. Furthermore, the 
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absorption is exceptionally well isolated from other spectral fea- 
tures. 

Spectroscopic observations in the near infrared (H/K band) 
would be a very valuable diagnostics. We hope that such observa- 
tions, taken simultaneously with spectra in th e optical, will becom e 
available for different SN Ib/Ic subtypes (e.g. iNomoto et al.l200ot) . 



6 SUMMARY AND CONCLUSIONS 

We used an NLTE spectral synthesis code to determine the ef- 
fect of different H and He masses on the spectra of stripped core- 
collapse SNe with low/moderate explosion energies and ejecta 
masses (10 51 erg, 1-3Mq). To this end, we first created best-fit 
models for SNe 2008ax (Type lib) and 19941 (Type Ic) These 
are in line with earlier models for the two SNe JSauer et alj[2006l : 
iTsvetkov et"ai]|2009l ; rMaurer et ai1l20ld) . We then calculated a se- 
quence of transitional models for which the envelope configura- 
tion is intermediate between the two extreme cases. From the se- 
quence models, we derived limits on the H and He content of low- 
mass/low-energy Type lb and Type Ic SNe, respectively. 

A small H mass (0.025 -O.O33A/0) is sufficient to make a 
strong Ha absorption show up. We speculate that at least some SNe 
lb show an Ha absorption component, blended with Si II A6355. 
This makes the spectral distinction between SNe lib and SNe lb 
somewhat fuzzy. Despite this, and the dependence of the H state 
on details in the models, we can however say that SNe lb with an 
ejecta mass of 2-3Mq should have an envelope containing less 
than 0.1 M Q of H. 

The maximum He mass which can be hidden in SNe Ic ac- 
cording to our models is ~0.06-0. 1471/q. When He lines appear, 
He I 2/j, absorption is an excellent diagnostic of He content. More 
SN Ib/Ic data in this spectral range are clearly needed. 

An assessment of the optical He I lines at AA6678, 7065 has 
however been sufficient to spot a difficulty wi th SN Ic progenitor 
models: Recent s tellar evolution calculations dGeorgy et alj|2009l ; 
lYoonetalj|2010l) generally predict a quite large He mass in SN 
Ilb/Ib/Ic progenitors (~0.3 -O.6M0, with few exceptions). "Ex- 
otic" models may therefore be needed to e xplain SN 19941-like , 
low-mass SNe Ic. Ilwamoto et al.l l fl994l) and lNomoto et al.l i ll 9941) 
proposed for this SN a binary scenario in which the companion is 
a neutron star in the last mass transfer phase. The rates of this pro- 
genitor channel may not be all too large; whether this is a problem 
or not depends on the exact rate of 19941-like SNe. 

It has to be verified whether more He can be hidden in the 
ejecta of more massive or strongly asymmetric SNe Ic. We there- 
fore plan to carry out similar studies on spectra of core-collapse su- 
pernovae which are more massive and more energetic or markedly 
aspherical. These studies will allow us to further explore how much 
H or He can be hidden with different configurations of the envelope 
(mixing, density profile). Together with the present work, we will 
obtain a comprehensive picture of the outer layers of SNe lib, lb 
and Ic, and of possible progenitor systems. 



ACKNOWLEDGEMENTS 

SH thanks N. Przybilla for very helpful discussions about atomic 
models and for a sample He I model, and I. Maurer for several 
fruitful discussions and routines for treating H atomic data. We are 
indebted to P. Storey for providing us with photoionisation rates 



for He I. ST acknowledges support by the Transregional Collabo- 
rative Research Centre TRR33 "The Dark Universe" of the German 
Research Foundation (DFG). We have made use of the CHIANTI 
(v6.0.1), TMAD (as available in 2010) and NIST (v3.1.5) databases 
for atomic data. CHIANTI (http://www.chiantidatabase.orgl is a 
collaborative project involving the NRL (USA), RAL (UK), MSSL 
(UK), the Universities of Florence (Italy) and Cambridge (UK), 
and George Mason University (USA). TMAD is available in the 
context of the German Astrophysical Virtual Observatory (GAVO) 
at |http://astro.uni-tuebingen.de/^rauc h/TMAD/TMAD.htmi] The 
NIST atomic spectra database is provided by the National Insti- 
tute of Science and Technology, Gaithersburg, MD, USA. It is ac- 
cessible (together with individual references for the atomic data) at 
http://www.nist.gov/pml/data/asd.cfm 



References 

Abbott D. C, Lucy L. B., 1985, ApJ, 288, 679 
Arcavi I., et al., 2011, ApJ, accepted, arXiv: 1106.3551 [astro- 
ph.CO] 

Arnett W. D., 1974, ApJ, 193, 169 

Axelrod T. S., 1980, PhD thesis, University of California, Santa 
Cruz 

Bashkin S., Stoner J. O., 1975. North-Holland Publ. Co., Amster- 
dam 

Branch D., Jeffery D. J., Young T. R., Baron E., 2006, PASP, 118, 
791 

Cann N. M., Thakkar A. J., 2002, J. Phys. B, 35, 421 
Cappellaro E., Mazzali P. A., Benetti S., Danziger I. J., Turatto 

M., della Valle M., Patat E, 1997, A&A, 328, 203 
Chornock R., et al., 2010, ArXiv: 1001.2775 
Chugai N. N., 1987a, Astrofizika, 26, 89 
Chugai N. N., 1987b, SvA, 13, L282 

Clocchiatti A., Wheeler J. C, Brotherton M. S., Cochran A. L., 

Wills D., Barker E. S., 1996, ApJ, 462, 462 
Dere K. P., Landi E., Mason H. E., Monsignori Fossi B. C, Young 

P. R., 1997, A&AS, 125, 149 
Dere K. P., Landi E., Young P. R., Del Zanna G., Landini M., 

Mason H. E., 2009, A&A, 498, 915 
Deuflhard P., 2004, Newton Methods for Nonlinear Problems: 

Affine Invariance and Adaptive Algorithms. Springer, Berlin 
Drake G. W. E, Morton D. C, 2007, ApJS, 170, 251 
Fernley J. A., Seaton M. J., Taylor K. T, 1987, J. Phys. B, 20, 

6457 

Filippenko A. V., et al., 1995, ApJ, 450, LI 1 + 

Georgy C, Meynet G., Walder R., Folini D., Maeder A., 2009, 

A&A, 502, 611 
Gould R. J., 1972, Physica, 60, 145 
Graham J. R., 1988, ApJ, 335, L53 
Habing H. J., Goldsmith D. W., 1971, ApJ, 166, 525 
Hachinger S., 201 1, PhD thesis, Technische Universitat Miinchen, 

Munich, Germany 
Hamuy M., et al., 2002, AJ, 124, 417 
Harkness R. P., et al., 1987, ApJ, 317, 355 
Hummer D. G., Rybicki G. B., 1985, ApJ, 293, 258 
Hummer D. G., Storey P. J., 1998, MNRAS, 297, 1073 
Iwamoto K, Nomoto K, Hoflich P., Yamaoka H., Kumagai S., 

Shigeyama T, 1994, ApJ, 437, LI 15 
Janev R. K, Smith J. J., 1993, Atomic and plasma-material inter- 
action data for fusion (Supplement to the journal Nucl. E), Vol. 
4. International atomic energy agency, Vienna 



© 2012 The Authors. Journal compilation © 2012 RAS, MNRAS 



He in SNe Ic, and H in SNe lb? 17 



Janka H., Langanke K., Marek A., Martfnez-Pinedo G., Miiller B., 
2007, Phys. Rept., 442, 38 

Khokhlov A., Mueller E., Hoeflich P., 1993, A&A, 270, 223 

Kozma C, Fransson C, 1992, ApJ, 390, 602 

Lach G., Pachucki K., 2001, Phys. Rev. A, 64, 042510 

Lamers H. J. G. L. M., Cassinelli J. P., 1999, Introduction to Stel- 
lar Winds. Cambridge University Press, Cambridge 

Leonard D. C, et al., 2006, Nature, 440, 505 

Leonard D. C, Filippenko A. V., 2005, in Turatto M., Benetti 
S., Zampieri L., Shea W., eds, AIP Conference Series Vol. 342, 
Spectropolarimetry of Core-Collapse SNe. AIP, Melville, p. 330 

Lotz W., 1967, ApJS, 14, 207 

Lucy L. B., 1991, ApJ, 383, 308 

Lucy L. B., 1999, A&A, 345, 211 

Lucy L. B., 2002, A&A, 384, 725 

MacFadyen A. I., Woosley S. E., 1999, ApJ, 524, 262 

Matzner C. D., McKee C. F, 1999, ApJ, 510, 379 

Maurer I., 2010, PhD thesis, Technische Universitat Miinchen, 
Munich, Germany 

Maurer I., Mazzali P. A., Taubenberger S., Hachinger S., 2010, 
MNRAS, 409, 1441 

Mazzali P. A., 2000, A&A, 363, 705 

Mazzali P. A., et al„ 2008, Sci, 321, 1185 

Mazzali P. A., Lucy L. B., 1993, A&A, 279, 447 

Mazzali P. A., Lucy L. B., 1998, MNRAS, 295, 428 

Mazzali P. A., Nomoto K., Cappellaro E., Nakamura T., Umeda 
H., Iwamoto K., 2001, ApJ, 547, 988 

Mezzacappa A., 2005, Annual Review of Nuclear and Particle 
Science, 55, 467 

Mihalas D., 1978, Stellar atmospheres, 2nd edn. W. H. Freeman 
and Co., San Francisco 

Mihalas D., Stone M. E., 1968, ApJ, 151, 293 

Nomoto K., Tominaga N, Tanaka M., Maeda K., Suzuki T., Deng 
J. S., Mazzali P. A., 2006, Nuovo Cimento B Serie, 121, 1207 

Nomoto K., Yamaoka H., Pols O. R., van den Heuvel E. P. J., 
Iwamoto K., Kumagai S., Shigeyama T., 1994, Nature, 371, 227 

Nowak U., Weimann L., 1991, Technical Report TR-91-10, A 
Family of Newton Codes for Systems of Highly Nonlinear Equa- 
tions. Konrad-Zuse-Zentrum fuer Informationstechnik Berlin 

Opal C. B., Beaty E. C, Peterson W. K., 1972, Atomic Data, 4, 
209 

Pastorello A., et al., 2008, MNRAS, 389, 955 
Piran T., 2005, Rev. Mod. Phys., 76, 1 143 

Ralchenko Y, Janev R. K., Kato T., Fursa D. V., Bray I., de Heer 
F. J., 2008, ADNDT, 94, 603 

Rauch T., Deetjen J. L., 2003, in Hubeny I., Mihalas D., Werner 
K., eds. Stellar Atmosphere Modeling, Vol. 288 of ASP Con- 
ference Series, Handling of Atomic Data. ASP, San Francisco, 
p. 103 

Richmond M. W., et al., 1996, AJ, 1 1 1, 327 

Roming P. W. A., et al., 2009, ApJ, 704, LI 18 

Sasaki M., Kosugi G., Ishigaki T., Maemura H., Aoki K., Ohtani 
H., 1994, PASJ, 46, L187 

Sauer D. N, Mazzali P. A., Deng J., Valenti S., Nomoto K., Filip- 
penko A. V, 2006, MNRAS, 369, 1939 

Schlegel D. J., Finkbeiner D. P., Davis M., 1998, ApJ, 500, 525 

Shigeyama T., Nomoto K., Tsujimoto T., Hashimoto M., 1990, 
ApJ, 361, L23 

Shigeyama T., Suzuki T., Kumagai S., Nomoto K., Saio H., Ya- 
maoka H., 1994, ApJ, 420, 341 
Smartt S. J., 2009, ARA&A, 47, 63 



Stehle M., Mazzali P. A., Benetti S., Hillebrandt W„ 2005, 

MNRAS, 360, 1231 
Swartz D. A., 1994, ApJ, 428, 267 
Taubenberger S., et al., 201 1, MNRAS, 413, 2140 
Tsvetkov D. Y, Volkov I. M., Baklanov P., Blinnikov S., Tuchin 

O., 2009, Peremennye Zvezdy, 29, 2 
Utrobin V. P., 1996, A&A, 306, 219 
Valenti S., et al., 201 1, MNRAS, p. 1220 
Woosley S. E., 1993, ApJ, 405, 273 
Xu Y, McCray R., 1991, ApJ, 375, 190 

Yokoo T, Arimoto J., Matsumoto K., Takahashi A., Sadakane K., 

1994, PASJ, 46, L191 
Yoon S., Woosley S. E„ Langer N, 2010, ApJ, 725, 940 



APPENDIX A: MC ESTIMATOR FOR THE 
FREQUENCY-DEPENDENT INTENSITY 

Here, we briefly comment on the implementation of the MC radi- 
ation field estimator ,l v mentioned in Sec. 12. II In order to evaluate 
the estimator in our simulations, packets are registered at each cell 
midpoint. The packet flux is binned on a sufficiently fine frequency 
grid; the result is then converted to a co-moving frame intensity and 
finally smoothed somewhat. At A < 1000 A, the recorded spectrum 
is relatively noisy. Therefore, we impose a continuous transition 
to a constant mean intensity of 10 _30 ergHz _1 s _1 sr~ 1 cm -2 to- 
wards smaller wavelengths. For A < 300 A, this normally exceeds 
the black-body intensity at the local electron temperature T e , and 
we let the intensity decline following B v (T e ). The exact values 
[10" 30 ergHz" 1 s~ 1 sr" 1 cm~ 2 ,B^(r e )] turned out to be irrelevant 
for the present work, as in our simulations most ionisations of He 
occur by nonthermal collision processes, and H is mostly ionised 
by photons with > 900 A (also from excited states), so that the 
assumed intensity does not control the rate. 



APPENDIX B: GAMMA-RAY DEPOSITION 
CALCULATIONS 

In our NLTE calculations, we take into account that atoms can be 
excited and ionised by Compton electrons, which are present owing 
to the gamma radiation from 56 Ni and 56 Co decay. In this appendix, 
we show how we determine what fraction D k of the local heating 
power He is deposited in different channels k, corresponding to 
atomic processes or to Coulomb scattering with ther mal electrons . 
Our implementation follows the principles devised by |Lucvl ( ll99 ll ). 

As a basic assumption, we consider the number of atoms in 
excited states to be small and insignificant, so that only excit ations 
and io nis ations from the gr ound state must be considered (cf. lLucvl 
Il99ll and lHachingerll201ll) . The deposition fraction Du is a func- 
tion of the initial energy of the fast electrons E. Assuming that all 
electrons which contribute to He have the same initial energy E, 
we can write the energy flow into channel k as: 

H k = D k {E) H c . 

For E ranging fr om ~ 1 keV up to some Me V, D k (E) becomes 
nearly constant dLucvl 1 199 it lAxelrod Il980h . Thus, for typical 
Compton electrons in a SN, we can write 

Hk = D+Hc, 

with an energy-independent value D~£ . 

We now describe how we compute the values D k (E) or . 
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The deposition process involves the creation of fast secondary elec- 
trons (mainly in ionisations). Therefore, a cascade of slower and 
slower electrons is produced. This is taken into a ccount solving 
an integrated vers ion of the Spencer-Fano equation dXu & McCravl 
ll99lULucvl[T99ll) . which gives the deposition fraction D k (E) for 
a process k at a source energy E in terms of the deposition fraction 
Dk(E') at all lower energies E' < E: 

ED k (E) = W kPk (E) + J2 Pj(E)-(E-W J )-D k (E-W j ) 

{j\E-W J ^E k , 
one e — -proc.} 

(E-Wt)/2 

+ Pi (E) ■ I dE s U l {E s )-\ i E s -D k (E 3 ) + 

{i\E-W^E k , { 
two e — -proc.} 

(E-E s -W l )-D k {E-E s -W^ (Bl) 

Here, i, j and k are indices of processes (channels), and p k (E) is 
the probability for process k to occur for the electron at energy E, 
which is proportional to the respective cross section. The energy 
loss (to atoms or thermal electrons) in process i\j\k is denoted as 
Wj\j\ k . Scatt ering on thermal electrons [for stopping numbers, see 
iGouldl l l 19721) 1 is, for simplicity, treated as a discretised process in 
which a nonthermal electron loses 5% of its energy to the plasma. 
This is not critical if the electrons ar e coming from large enough en- 
ergies faabing & Goldsmith|[T97lh . The different terms contribut- 
ing to ED k (E) are as follows: The W k p k (E) term accounts for the 
direct deposition into the Channel k. The other two terms account 
for contributions from electrons which first trigger a process j or 
i and then deposit energy into channel k. Because of such contri- 
butions, the deposition fraction at some energy E is related to that 
at smaller energies. The latter of these terms specifically takes into 
account ionisation processes which produce secondary electrons. 
Hi(E s ) is the probability distribution for the energy of the sec- 
ondary electron E a , and integrates up to one. Following lOpal et"al] 
dl972l) . we assume Ui(E s ) oc 1/(E 2 S + W 2 ), where W is 13.8 eV 
for He I and corresponds to the ionisation energy for other species. 

Since only atomic processes from the ground state are taken 
into account, non-thermal electrons with an energy below some 
threshold can only scatter on thermal electrons. With initial values 
for D k set up accordingly, we integrate the system l lB 1 1 up to an en- 
ergy of 10 keV. The deposition fractions then obtained are adopted 
as . We have tabulated these values for different H/He mixes, 
ionisation states, fractions nm+Ht\/n e (where all atoms are as- 
sumed to be in the ground state) and absolute va lues of n e (although 
the d ependence on the latter is very weak, cf. iKozma & Franssonl 
1 1992b . The tables are read in by the NLTE module depending on 
the actual local composition. If elements other than H and He are 
present, we first approximate the rates assuming an ?iH+He as large 
as the total number density of atoms n tot . Then, we scale down the 
value for every process by a factor: 

p _ Sj = H,He ZjTli 
all species j J J 

where Zi are atomic numbers and rii are number densities. This 
choice of the scaling factor corresponds to the crude assumption 
that the energy deposition by fast electrons into a certain chemical 
species is proportional to the number density of that species and 
to the number of electrons in one atom. In section [5T1 we show 
that our results - at least for He - are surprisingly insensitive to 
the exact exci tation rates, so that this approximation [for "exact" 
values, cf. e.g. lLotzl dl967h l is sufficient for the present study. 



APPENDIX C: INPUT PARAMETERS FOR THE 
SPECTRAL SYNTHESIS CODE 

In Table ICll we give the code input parameters for our SN 2008ax 
and SN 19941 models. The table with its annotations illustrates the 
one-to-one shell correspondence we have established between the 
abundance shells in the two models for the SNe (cf. Sec. 14.1.21 . 
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Table CI. Abundance shells and code input parameters for the models for SN 2008ax and SN 19941. The models are specified such that the abundance shells 
in SN 2008ax and SN 19941 correspond to one another (i.e. each column in the upper part of the table corresponds to the same column in the lower part), ^ibord 
is the lower border for each shell. Shells with epochs given correspond to observed spectra (i.e. their i% OI i is the photospheric velocity at that time). Shells 
with epochs "o" or "oo" are outer shells set up to match the t = 1 6.0 d spectrum. Shells with epoch "-" are inserted in order to enable a one-to-one mapping 
of SN 2008ax to SN 19941 (see also footnotes to the table). We emphasize that the abundance structure in each of the two models, as far as it contributes to the 
formation of the respective spectra, is effectively that of a six-zone model (i.e. some of the eleven zones in each model have equal abundances or are irrelevant 
for the spectra). 

SN 2008ax 



epoch [d rel. to B max.] 










t [d rel. to explosion] 










1 










"lbord [km s ] 


390 


670 


1500 


2500 


X(H) [%] 


0.0000 


0.0000 


0.0000 


0.0000 


X(He) [%] 


0.0000 


0.0000 


0.0000 


1.0000 


X(C) [%] 


0.1500 


0.5000 


1.5000 


5.0000 


X(O) [%] 


47.1097 


48.6097 


56.8597 


63.2158 


X(Ne) [%] 


0.0000 


0.5000 


6.0000 


6.0000 


X(Na) [%] 


0.1000 


0.2500 


4.0000 


4.0000 


X(Mg) [%] 


1.5000 


2.5000 


12.0000 


10.0000 


X(Si) [%] 


8.0000 


11.0000 


10.0000 


5.0000 


X(S) [«] 


2.5000 


4.0000 


5.0000 


2.0000 


X(Ca) [%] 


0.4200 


0.4200 


0.4200 


0.2700 


X(Ti) [%] 


0.0200 


0.0200 


0.0200 


0.0800 


X(Cr) [%] 


0.2000 


0.2000 


0.2000 


0.2000 


X(Fe ) c [%] 


2.5000 


4.0000 


0.5000 


0.0333 


X( 56 Ni ) c [%] 


37.5000 


28.0000 


3.5000 


3.2000 



+ 21.7 




+ 11.7 


+ 3.2 


- 2.9 





oo 


40.60 




30.60 


22.05 


16.00 






8.390 




8.530 


8.745 


8.695 






3100 


3300 


6000 


7300 


8400 


10800 


13000 


0.0000 


0.0000 


0.0000 


9.0000 


15.0000 


20.0000 


40.0000 


45.0000 


45.0000 


51.0000 


57.0000 


60.0000 


69.5000 


56.0000 


17.0000 


17.0000 


15.0000 


10.0000 


8.0000 


5.0000 


1.5000 


17.1917 


17.1917 


14.4167 


10.6717 


6.2650 


2.3417 


0.9791 


3.0000 


3.0000 


2.5000 


2.0000 


1.5000 


0.8000 


0.4000 


1.5000 


1.5000 


1.5000 


1.0000 


1.0000 


0.6000 


0.3000 


6.0000 


6.0000 


6.0000 


4.0000 


3.0000 


0.6000 


0.3000 


5.0000 


5.0000 


5.0000 


3.0000 


2.5000 


0.7500 


0.3750 


2.0000 


2.0000 


2.0000 


1.2500 


1.0000 


0.3000 


0.1000 


0.0500 


0.0500 


0.0250 


0.0250 


0.0150 


0.0060 


0.0003 


0.0050 


0.0050 


0.0050 


0.0050 


0.0050 


0.0003 


0.0003 


0.0200 


0.0200 


0.0200 


0.0150 


0.0150 


0.0020 


0.0020 


0.0333 


0.0333 


0.0333 


0.0333 


0.0333 


0.0333 


0.0333 


3.2000 


3.2000 


2.5000 


2.0000 


1.6667 


0.0667 


0.0100 



SN 19941 



epoch [d rel. to B max.] + 28.6 

£ [d rel. to explosion] 40.60 

Ig(^gl) 7.912 

"lbord [kms -1 ] 910 

X(H) [%] 0.0000 

X(He) [%] 0.0000 

X(C) [%] 0.1500 

X(O) [%] 0.1097 

X(Ne) [%] 0.0000 

X(Na) [%] 0.1000 

X(Mg) [%] 1.5000 

X(Si) [%] 8.0000 

X(S) [%] 2.5000 

X(Ca) [%] 0.4200 

X(Ti) [%] 0.0200 

X(Cr) [%] 0.2000 

X(Fe ) c [%] 7.0000 

X( 56 Ni ) c [%] 80.0000 



+ 18.6 


+ 10.1 


+ 4.0 








_b 




b 


b 


b 


30.60 


22.05 


16.00 
















8.122 


8.430 


8.770 
















1400 


3100 


8900 


15000 


18500 












0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


9.0000 


15.0000 


20.0000 


40.0000 


0.0000 


0.0000 


1.0000 


5.0000 


5.0000 


65.7058 


66.4508 


67.3575 


70.5309 


56.4002 


0.5000 


1.5000 


5.0000 


38.0000 


44.0000 


15.0000 


10.0000 


8.0000 


5.0000 


1.5000 


5.6097 


43.8597 


70.6658 


32.8059 


32.1225 


14.9167 


10.6717 


6.2650 


2.3417 


0.9791 


0.5000 


6.0000 


6.0000 


20.0000 


17.0000 


2.5000 


2.0000 


1.5000 


0.8000 


0.4000 


0.2500 


4.0000 


4.0000 


2.0000 


1.0000 


1.0000 


1.0000 


1.0000 


0.6000 


0.3000 


2.5000 


12.0000 


10.0000 


2.0000 


0.7500 


0.7500 


0.7500 


0.7500 


0.6000 


0.3000 


11.0000 


10.0000 


0.7000 


0.0500 


0.0500 


0.0500 


0.0500 


0.0500 


0.0500 


0.0500 


4.0000 


5.0000 


0.2500 


0.0250 


0.0250 


0.0250 


0.0250 


0.0250 


0.0250 


0.0250 


0.4200 


0.4200 


0.2700 


0.0167 


0.0001 


0.0001 


0.0001 


0.0001 


0.0001 


0.0001 


0.0200 


0.0200 


0.0800 


0.0004 


0.0004 


0.0004 


0.0004 


0.0004 


0.0003 


0.0003 


0.2000 


0.2000 


0.2000 


0.0020 


0.0020 


0.0020 


0.0020 


0.0020 


0.0020 


0.0020 


5.0000 


2.0000 


0.0333 


0.0333 


0.0333 


0.0333 


0.0333 


0.0333 


0.0333 


0.0333 


70.0000 


15.0000 


1.8000 


0.0667 


0.0167 


0.0167 


0.0167 


0.0167 


0.0167 


0.0100 



a In our model SN 2008ax, we had to introduce zones below the i — 40.6 d photosphere. The abundances have generally been chosen 
as in the corresponding zones in SN 19941. Some differences appear, as we had to adapt the Fe-group abundances because the light curve 
of SN 2008ax required a shallower 56 Ni distribution than what we find in SN 19941. 

b In the SN 19941 model we formally added five shells at mass coordinates "outside" the density model (i.e. at M > 0.9Mq). These 
shells do not have any impact on the SN 19941 calculations. They just correspond ( in mass coordinates) to the outer shells in SN 2008ax 
and inherit the respective abundances with a few exceptions: the abundances of IME and Fe-group elements have been limited to those in 
the outermost shells of our original SN 19941 model (o, oo). 

c The abundances of Fe, Co and Ni in our models are assumed to be the sum of j6 Ni and the elements produced in the decay chain 
( 56 Co and 56 Fe) on the one hand, and directly synthesised / progenitor Fe on the other hand. Thus, they are conveniently given in terms 
of the 56 Ni mass fraction at t — [X ( 56 Ni)o], the Fe abundance at t — [X (Fe)oL and the time from explosion onset t. 



© 2012 The Authors. Journal compilation © 2012 RAS, MNRAS 



